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We simulate the central reactions of nearly symmetric, and asymmetric systems, for 
the energies at which the maximum production of IMFs occurs (E^ fc ).This study is 
carried out by using hard EOS along with cugnon cross section and employing MSTB 
method for clusterization. We study the various properties of fragments. The stability 



of fragments is checked through persistence coefficient and gain term. The information 
about the thermalization and stopping in heavy-ion collisions is obtained via relative 
momentum, anisotropy ratio, and rapidity distribution. We find that for a complete 
stopping of incoming nuclei very heavy systems are required. The mass dependence of 
various quantities (such as average and maximum central density, collision dynamics as 
well as the time zone for hot and dense nuclear matter) is also presented. In all cases 
(i.e., average and maximum central density, collision dynamics as well as the time zone 
for hot and dense nuclear matter) a power law dependence is obtained. 

1 Introduction 

The heavy-ion collisions at intermediate energies are excellent tool to study the nuclear 
matter at high density and temperature. At high excitation energies, the colliding nuclei 
compress each other as well as heat the matter [1-3]. This leads to the destruction of 
initial correlations, which in turn makes the matter homogeneous and one can have global 
stopping. The global stopping is defined as the randomization of one-body momentum 
space or memory loss of the incoming momentum. The degree of stopping, however, may 
vary drastically with incident energies, mass of colliding nuclei and colliding geometry. 
The degree of stopping has also been linked with the thermalization (equilibrium) in 
heavy- ion collisions. More the initial memory of nucleons is lost, better it is stopped. 

The fragmentation of colliding nuclei into several pieces of different sizes is a complex 
phenomenon. This may be due to the interplay of correlations and fluctuations emerging 
in a collision. Several studies, in literature, have been made to check the fragmentation 
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pattern. Fragmentation pattern has been reported to depend on the size of the colliding 
nuclei, incident energy as well as impact parameter [1,4-6]. Dhawan et al. [7] studied 
the degree of stopping reached in intermediate energy heavy-ion collisions. They found 
that degree of stopping decreases with increase in impact parameter as well as at very 
high energies. They suggested that the light charged particles (LCPs) (2<A<4), can 
be used as a barometer for studying the stopping in heavy-ion collisions. As lighter 
fragments mostly originate from midrapidity region whereas intermediate mass fragments 
(IMFs) originate from surface of colliding nuclei and can be viewed as remnants of the 
spectator matter. On the other hand, Sood and Puri [8] studied the thermalization 
achieved in heavy-ion collisions in terms of participant-spectator matter. They found that 
participant-spectator matter depends crucially on the collision dynamics as well as history 
of the nucleons and important changes in the momentum space occur due to the binary 
nucleon-nucleon collisions experienced during the high dense phase. The collisions push 
the colliding nucleons into midrapidity region responsible for the formation of participant 
matter. This ultimately leads to thermalization in heavy-ion collisions. Vermani et al. [8] 
used rapidity distribution of nucleons to characterize the stopping and thermalization of 
the nuclear matter. They found that nearly full stopping is achieved in heavier systems 
like 197 Au+ 197 Au whereas in lighter systems a larger fraction of particles is concentrated 
near target and projectile rapidities, resulting in a broad Gaussian shape. The lighter 
systems, therefore, exhibit larger transparency effect, i.e., less stopping. Puri et al. [9] 
studied the non-equilibrium effects and thermal properties of heavy-ion collisions. They 
found that the heavier masses are found to be equilibrated more than the lighter systems. 
Recently, Sisan et al. [10] studied the emission of IMFs from central collisions of nearly 
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symmetric systems using a 47r-array set up, where they found that the multiplicity of IMFs 
shows a rise and fall with increase in the beam energy. They observed that E^f (the 
energy at which the maximum production of IMFs occurs) increases linearly with the 
system mass, whereas a power-law (oc A r ) dependence was reported for peak multiplicity 
of IMFs with power factor r = 0.7. Though percolation calculations reported in that paper 
failed to explain the data, subsequent calculations using QMD model [11] successfully 
reproduced the data over entire mass. One is, therefore, interested to understands how 
nuclear dynamics behaves at this peak energy i.e. whether linear increase reported for 
the multiplicity of fragments remains valid for other observables or not. We here plan 
to investigate the degree of stopping reached and other related phenomena in heavy- 
ion reactions at peak center-of-mass energies. We also check system size dependence of 
average and maximum central density, collision dynamics as well as the time zone for hot 
and dense nuclear matter at peak center-of-mass energy. 

This study is made within the framework of the quantum molecular dynamics model, 
which is described in detail in Refs. [1,4, 11-18]. 

2 Results and Discussion 

For the present study, we simulate the central reactions (b =0.0 fm) of 20 Ne+ 20 Ne, 
40 Ar+ 45 Sc, 58 Ni+ 58 Ni, 86 Kr+ 93 Nb, 129 Xe+ 118 Sn, 86 Kr+ 197 Au and 197 Au+ 197 Au at the in- 
cident energies at which the maximal production of intermediate mass fragments (IMFs) 
occurs. These read approximately 24, 46, 69, 77, 96, 124, and 104 AMeV, respectively for 
the above mentioned systems [11]. Note that at lower incident energies phenomena like 
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fusion, fission and cluster decay are dominant [19]. Here we use hard (labeled as Hard) 
equation of state along with energy dependent Cugnon cross section (o"{^ ee ) [15]. The 
reactions are followed till 200 fm/c. The phase-space is clusterized using minimum span- 
ning tree method with binding energy check (MSTB). The MSTB method is an improved 
version of normal MST method [16]. Firstly, the simulated phase-space is analyzed with 
MST method and pre-clusters are sorted out. Each of the pre-clusters is then subjected 
to binding energy check [11, 16]: 

1 Nf (r> _p c -"M2 i Nf 
&=Nf B — 2m Nf + 2 £ V *<*> r ^ < ^ 

We take E bind = -4.0 MeV if > 3 and E bind = 0.0 otherwise. Here N f is the number 
of nucleons in a fragment and P^/' is center-of-mass momentum of the fragment. This is 
known as Minimum Spanning Tree method with Binding energy check (MSTB) [11, 16]. 
Note that nucleons belong to a fragment if inequality (1) is satisfied. The fragments formed 
with the MSTB method are more reliable and stable at early stages of the reactions. 

One of the important aspects in fragmentation is the stability of fragments as well 
as surrounding nucleons of a fragment. The change in the nucleon content of fragments 
between two successive time steps can be quantified with the help of persistence coefficient 
[17,18]. 

Let the number of pairs of nucleons in cluster C at time t is xc(t)=0.5*N c (N c -l)- At 
the time At later, it is possible that some of the nucleons belonging to cluster C have left 
the cluster and are the part of another cluster or are set free or others may have entered 
the cluster. Now, let N Cd be the number of nucleons that have been in the cluster C at 
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time t and are at t + At in cluster D. We define 



$c(t + At) = J2 0-5 * N Cd (N Cd - 1). 



(2) 



D 



Here the sum runs over all fragments D present at time t + At. 



The persistence coefficient of cluster C can be defined as [17, 18]: 



Pc(t + —) = ^c(t + At)/ X c(t). 



(3) 



The persistence coefficient averaged over an ensemble of fragments is defined as: 



where N/ t is the number of fragments present at time t in a single simulation. The quantity 
is then averaged over a large number of QMD simulations. The stability of a fragment 
between two consecutive time steps can be measured through persistence coefficient. If 
fragment does not emit a nucleon between two time steps, the persistence coefficient is 
one. On the other hand, if fragment disintegrates completely, the persistence coefficient 
will be zero. If we remove one nucleon from a fragment C, the persistence coefficient is 
P c (t + At/2) = (N c -2)/N c i.e., 0.333 for N c = 3 and 0.8 for N c = 10. For example for 
mass 10, when one nucleon is emitted we have two entities at later time step consisting 
of free nucleon and fragment with mass 9. The Pc(t + At/2) is the contribution from 
all such entities existing at later times. It, then, measures the tendency of the members 
of given cluster to remain together. In fig. 1, we display the persistence coefficient 
for various fragments i.e., light charged particles (LCPs) (2<A<4), intermediate mass 
fragments (IMFs) (5<A<44) as well as heavy mass fragments (HMFs) (10<A<44). The 
various lines have been defined in the caption. It is clear from fig. that the saturation 
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value of persistence coefficient is slightly higher in case of LCPs as compared to heavier 
fragments. One can conclude that the final fragments are formed after 130 fm/c when 
this coefficient is 0.8. Before that time there is a strong exchange of nucleons between 
the fragments. The number of medium and intermediate size fragments increases because 
the largest fragment falls finally into this mass bracket. The persistent coefficient reaches 
its asymptotic value later due to the interaction between fragments as well as between 
fragments and free nucleons. Due to this interaction nucleons are sometimes absorbed or 
emitted from the fragments. This process changes the details but not the general structure 
of the fragmentation pattern. 

The persistence coefficient tells about the stability of different fragments between two 
successive time steps. But it does not provide any information whether a fragment has 
swallowed some nucleons or not. To check this, we use a quantity called "Gain" [18]. The 
Gain represents the percentage of nucleons that a fragment has swallowed between two 
consecutive time steps. Let N£ be the number of nucleons belong to a fragment a at time 
t. Let be the number of nucleons which were in cluster a at time t and are in cluster 
f5 at time t + At. The Gain is defined as: 



r] = 0.0, 0.5, and 1.0 if < 0.5 Nj, = 0.5 Nj and > 0.5 Nj, respectively. 
Naturally, a true Gain for a fragment a is only if its nucleons constitute at least half of 
the mass of new fragment (3. The Gain term will tell us whether the interactions among 
fragments have ceased to exist or not. 

In fig. 2, we display the gain term for LCPs, HMFs, and IMFs. As discussed earlier, 
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the value of persistence coefficient is slightly higher in case of LCPs. Therefore, gain term 
will be smaller for LCPs as shown in fig. 2. As evident from fig., in case of heavier systems 
the gain term has higher value because of the large nucleon-nucleon interactions. 

The quantities which are closely related to the degree of thermalization are relative 
momentum (K R ) and anisotropy ratio (R a ). The average relative momentum of two 
colliding fermi spheres is defined as [3,9,20]: 

(K R ) = {\P P (T,t)-P T (T,t)\/h), (6) 

where 

^i ^y . (7) 

Here Pj and pj are the momentum and density experienced by j th particle and k stands 
for either target or projectile and r refers to a space point in central sphere of 2 fm radius 
to which all calculations are made. The (K R ) is an indicator of local equilibrium because 
it depends on the local position r . 

The second quantity is anisotropy ratio which is defined as [3, 7, 9, 20]: 

(Ra) = V . (8) 

The anisotropy ratio (R a ) is an indicator of global equilibrium of the system because 
it represents the equilibrium of the whole system and does not depend upon the local 
positions. The full global equilibrium averaged over large number of events will correspond 
to (R a ) = 1. 

In figs. 3a and 3b, we display, respectively, (K R ) and (R a ) ratio as a function of time 
for different system masses. The initial value of relative momentum increases whereas of 
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the anisotropy ratio decreases with mass of the system since E^ fc increases with increase 
in the system mass. It is interesting to see that the relative momentum is large at the start 
of the reaction, and finally at the end of the reaction, the value of (K R ) is nearly zero. This 
means that at the end of the reaction, the local equilibrium is nearly reached. However, 
the saturation time is nearly the same throughout the mass range. It is clear from the 
fig. 3b anisotropy ratio changes to a greater extent during the high density phase. Once 
the high density phase is over, no more changes occur in thermalization. Interestingly, 
the heavier nuclei are able to equilibrate more than the lighter nuclei. This is because of 
the fact that the number of collisions per nucleon for the 197 Au+ 197 Au reaction is larger 
than for the 58 Ni+ 58 Ni reaction. 

The rapidity distribution is also assumed to give information about the degree of 
thermalization achieved in heavy-ion reactions. The rapidity distribution of i th particle is 
defined as [7]: 

W 2 E(i)-p,(i) 

Here E(i) and p z (i) are, respectively, the total energy and longitudinal momentum of i th 
particle. Naturally, for a complete equilibrium a single Gaussian shape peak is expected. 
In fig. 4, we display the rapidity distribution of free-nucleons, LCPs, HMFs as well as 
IMFs. Rapidity distribution of all types of fragments indicate that heavier systems are 
better thermalized as compared to lighter ones. For lighter nuclei, we get relatively flat 
distribution. The effect is more pronounced for different kinds of fragments as compared 
to free-nucleons. 

In fig. 5, we display the system size dependence of the maximal value of average 
density (p avg ) (solid circles) and maximum density (p ma:r ) (solid stars). Lines represent 
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the power law fitting (oc A T ). The maximal values of (p av9 ) and (p ma:r ) follow a power 
law (oc A T ) with r being 0.08±0.02 for the average density (p a '° 9 } and 0.034±0.008 for 
maximum density (p max ) i.e., a slight increase in density occurs with increase in the size 
of the system. This is because, E^ fc increases with the size of the system. 

In fig. 6, we display the time of maximal collision rate (open stars) and average density 
< Pavg > max (solid circles) as a function of the total mass of the system. Interestingly, 
both quantities show a nearly mass independent behavior which shows that E^ fc increases 
with the mass of the system in such a way that maximal collision rate and maximal density 
is achieved at the same time throughout the mass range. The power factor being equal 
to -0.017±0.038 for maximal time of collision rate and 0.015±0.07 for maximal time of 
average density. 

Apart from the maximal quantities, another interesting quantity is the dense zone at 
the peak energy. In fig. 7, we display the time interval for which p avg > p Q (solid circles) 
and p avg > p /2 (open stars). Again both quantities follow a power law behavior with r = 
0.15±0.05 and r— -0.04±0.06, respectively, for p avg > p Q and p avg > p /2. This indicates 
that the time duration for which p avg is greater than the normal nuclear matter density 
increases with the mass of the system. 

The system size dependence of the (allowed) nucleon-nucleon collisions (solid squares) 
is displayed in fig. 8. The results are displayed at 200 fm/c where the matter is di- 
luted and well separated. The nucleon-nucleon collisions increase with the system size. 
This enhancement can be parametrized with a power law proportional to A r with r = 
1.28±0.054. At fixed incident energy nucleon-nucleon collisions should scale as A. This 
has been tested by Sood and Puri [21]. Here power factor is greater than one since with 
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increase in mass of the system E£^ fc also increases. 

3 Summary 

In the present study, we have simulated the central reactions of nearly symmetric, and 
asymmetric systems, for the energies at which the maximum production of IMFs occurs 
(E^ fc ), using QMD model. This study is carried out by using hard EOS along with 
cugnon cross section and employing MSTB method for clusterization. We have studied 
the various properties of fragments. The stability of fragments is checked through persis- 
tence coefficient and gain term. The information about the thermalization and stopping 
in heavy-ion collisions is obtained via relative momentum, anisotropy ratio, and rapid- 
ity distribution. We found that for a complete stopping of incoming nuclei very heavy 
systems are required. The mass dependence of various quantities (such as average and 
maximum central density, collision dynamics as well as the time zone for hot and dense 
nuclear matter) is also presented. In all cases (i.e., average and maximum central density, 
collision dynamics as well as the time zone for hot and dense nuclear matter) a power law 
dependence is obtained. 
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Figure 1: The persistence coefficient as a function of time for LCPs, HMFs, and IMFs. 
Dashed, dotted, dash-dotted and solid lines are for 58 Ni+ 58 Ni, 86 Kr+ 93 Nb, 129 Xe+ 118 Sn 
and 197 Au+ 197 Au, respectively. 
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Figure 2: The Gain term as a function of time for LCPs, HMFs, and IMFs. Lines have 
same meaning as in fig. 1. 
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Figure 3: The time evolution of (a) relative momentum and (b) anisotropy ratio. Lines 
have same meaning as in fig. 1. 
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Figure 4: The rapidity distribution, dN/dY, as a function of reduced rapidity, Y c m /Y beam . 
Lines have same meaning as in fig. 1. 
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Figure 5: The maximal value of the average density (p av9 ) max (upper part) and maximum 
density {p max ) m ax (lower part) as a function of the composite mass of the system. The 
solid lines are the fits to the calculated resits using A r obtained with \ 2 minimization. 
The average is done over all space points on a sphere of 2 fm radius at center-of-mass. 
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Figure 6: The time of maximal value of collision rate (open stars) and average density 
(solid circles) as a function of composite mass of the system. The dashed and solid lines 
represent the x 2 fits with power law A r . 10 
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Figure 7: The time zone for p avg > p Q (upper part) and forp max > p Q (lower part) as a 
function of composite mass of the system. The solid lines represent the the \ 2 fits with 
power law A T . 9n 




Figure 8: The total number of allowed collisions versus composite mass of the system. 
The solid line represents the \ 2 fits with power law A r . 

21 



